Log-periodic Oscillations for Biased Diffusion on Random Lattice 
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ABSTRACT: Random walks with a fixed bias direction on randomly diluted 
cubic lattices far above the percolation threshold exhibit log-periodic oscillations in 
the effective exponent versus time. A scaling argument accounts for the numerical 
results in the limit of large biases and small dilution and shows the importance 
of the interplay of these two ingredients in the generation of the log-periodicity. 
These results show that log-periodicity is the dominant effect compared to previous 
predictions of and reports on anomalous diffusion. 

Diffusion on percolating lattices has a long history [1-4]. Here the sites of a large 
lattice are randomly initialized as being accessible, with probability p, or forbidden, 
with probability 1 — p', and random walker diffuse on the accessible sites only. In 
biased diffusion [5,6], with a probability B the random walk moves in one fixed 
direction, while with probability 1 — B it moves as usual to one randomly selected 
nearest neighbor; in both cases the move is allowed only if the neighbor is accessible. 
This fixed bias corresponds to an external electric field and is different from the bias 
with changing direction depending on the flow through the random network [7]. 

Numerical searches for a possible phase transition as a function of the bias B were 
hampered by strong variations of the effective exponent k in 

<r2>oct'^ (1) 

as a function of time t, for the mean square displacement r. In particular, Seifert 
and Suessenbach [8] found smooth variations, if k was plotted versus the logarithm 
of time. Log-periodic oscillations were predicted before for a special one-dimensional 
diffusion problem [9] while ref.8 investigated three dimensions. 

In the meantime, both computer technology and theoretical understanding have 
advanced appreciably. Log-periodic oscillations were seen in many cases, as reviewed 
recently [10], including the stock market crashes on Wall Street in 1929 and 1987 
[11] and in 1997 [12]. The physical mechanism for log-periodicity obtained in [9] is 
the exponential dependence of the trapping time on the length of the traps that are 
multiples of the lattice mesh size. The resulting intermittent random walk is thus 
punctuated by the successive encounters with larger and larger clusters of trapping 
sites. For strong bias B in three dimensions, we might expect a somewhat similar 
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mechanism as discussed below. We were thus motivated to repeat the work of ref.[8] 
for longer times, in search for periodicities like sin(log(t)). 

Each of the 512 processors of the Cray-T3E deals with about 3 million steps per 
second, roughly the same speed as the supercomputer CDC Cyber 205 used in ref.[8]. 
We could also simulate larger simple-cubic lattices (301 x 301 x 301) even though we 
wasted a whole 32-bit word for each lattice site (instead of one bit only) to simplify 
and accelerate the program (93 Fortran lines, available from the first author). We 
averaged over about 500 lattices, with 1000 walkers on each for short times and 3 for 
long times up to 4 x 10^°. 

If the bias equals B — 1, the motion is one-dimensional, and thus because of 
ref.[9] the chances for log-periodic oscillations are best for high B. Indeed, ref. [8] 
found clearer oscillations for B = 0.99 than for B = 0.9. Fig.l shows our results for 
B = 0.98 giving very clearly three maxima in k vs. log(t) near t = 10^, 2 10^, 4 10®. 
Fig. 2 shows our results for B — 0.999 giving two maxima in k vs. log(t) near 
t — 10^, 10^. Fig.3 shows our results for B — 0.95 giving many more maxima but 
with much smaller amplitude which arc thus harder to distinguish from noise. As a 
function of the bias 0.9 < B < 0.999 at fixed concentration p = 0.725 we found slower 
oscillations for larger B, while for smaller B they became more rapid. For smaller 
-B's, the oscillations become masked by the random fluctuations. We could probably 
attempt to measure them by a suitable averaging procedure [10,13] and leave this for 
the future. 

Table 1 : The table summarizes our estimates for the scaling factors A defined as 
the ratio of the times of the successive maxima. 
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We now present the scaling argument proposed to account for these results, that 
should be correct in the limit 5 — > 1 and 1 — p — > but non zero. We expect its 
domain of validity to be much broader than these limits would indicate, in a way 
similar to what we found for the scaling argument developed in [14] (see also [10]) to 
retrieve the exact analytical results derived in [9]. The mechanism is best understood 
by noticing that the biased random walk along x is the same as an unbiased random 
walk in the presence of a linearly decreasing potential 

V{x) = -bx , (1) 

where the slope b is such that the ratio of the forward over backward motion rate 
parallel to the bias is correctly evaluated. In the definition involving B, the forward 
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rate is + B = ^^^^ (the unbiased contribution plus the bias) and the backward 
rate is ^-g^- Equating their ratio to exp[/3 h] yields h in (|I|), where /3 is the inverse 
temperature in the proposed analogy : 

h (X B for small B and h oc log for B close to one. 

As the random walker makes its way through the system, it will eventually en- 
counter a local well, defined by a connected peninsula surrounded by empty sites 
(this can also be called a dead end) such that it has to retrace its steps backward to 
escape. Call n the depth of this well. The typical time of residence of the walker in 
this well (trapping time) is proportional to the Arrhenius factor exp[/3 h n\. This is 
the mechanism that converts a linear increase (in n) into an exponential increase in 
trapping times and thus produces log-periodic oscillations. Indeed, at short times, it 
is more probable that the smallest well be first encountered. This is due to the fact 
that such trapping wells occur with small probability (see below), denoted Pn- The 
smallest wells are of depth n = 1 and lead to a trapping time cx exp[/3 h ] = ^^^^ . 
The next most severe trapping wells are of depth n = 2 and lead to a typical trapping 
time oc exp[2(3 b ] = [ \^^^ ]^- And so on. We see a hierarchy of time scales [ ^l^jf ]" 
resulting from the succession of worst encountered trapping wells. Note that after it 
has encountered a well of depth n, the walker will in general encounter several wells 
of the same or smaller depths before being trapped in a well of depth n + 1. This 
does not destroy the argument because only the worst trapping well encountered until 
now dominates the waiting times and thus the log-periodic modulations. This is com- 
pletely similar to what was found in the one-dimensional case [9,10]. Note that we 
also predict that only a finite number rimax of oscillations will be observed in a finite 
system of linear size L. This simply stems from the finiteness of the deepest well rimax 
in a finite system, given by the condition Pn,nax L'^ ^ ^ (J- Machta and A. Aharony, 
priv.comm.). If F„ is exponentially decreasing with n, we see that Umax ~ logi^, 
while it is even smaller Umax ~ [logL]~ according to our suggestion below. This 
size effect has been verified qualitatively for B = 0.98, with the same concentration 
p = 0.725 and sizes L = 51, 21, 11 and 5. 

Our theory teaches us that only the depth of the wells control the leading behav- 
ior of the trapping time scales. Their transverse sizes bring in only subleading power 
law corrections to the exponential Arrhenius dependence. This is the fundamen- 
tal reason for the existence and robustness of the observed log-periodic oscillations, 
notwithstanding the many random configurations that the trapping wells can take. 
Notice that, as in all other known cases [10], the log-periodic oscillations rely on the 
existence of an ultraviolet (or infrared) cut-off, here the mesh size. 

We have reported in the last column of table 1 the value of the scaling factor 

_ tn+l _ 1 + 5B 

derived from (0). The agreement is quite good for the largest bias S — 1 for which the 
first smallest wells dominate completely the dynamics. For smaller S's, corrections 
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include the effect of all subleading time scales coming from all well shapes of a given 
depth. There is also a larger uncertainty in the determination of the scaling ratio A 
from the figures as the log-periodic oscillations are weaker and harder to distinguish 
from statistical noise. The small ratios seen for the smaller S's in the figures could 
also be higher harmonics. 

This is all we need to account for the observed log-periodicity. Let us now briefly 
comment on the distribution P„ of trapping well sizes of depth n. A trapping well 
of depth n has a typical transverse size n (isotropy). It can be constructed in a way 
similar to that leading to the distribution of finite cluster size above percolation. In 
this latter finite cluster s is constructed from the infinite incipient percolation 

network by introducing of the order of s~d~ empty perimeter sites. This leads to the 
stretched exponential distribution exp[— a s^~^^'^] known to be exact [15]. In order 
to construct a trapping well of depth n, one needs also to introduce of the order 
of s^~^^'^ empty perimeter sites. This would predict the distribution exp[— a' n'^"^], 
which decays faster than an exponential in dimensions larger than two. This is in 
contrast to [6] which used an exponential distribution to predict a transition from 
drift to no drift (vanishing velocity) above a critical value of the bias B. Our results 
show that this effect is very hard to detect numerically if true and that the main 
effect of the trapping wells is by far in creating a hierarchy of time scales t„ given by 
d^) leading to log-periodic oscillations that completely dominates the dynamics. 

To summarize, damped oscillations, periodic in the logarithm of the time, have 
been discovered for strongly biased diffusion on disordered three-dimensional lattices; 
these log-periodic oscillations stem from an Arrhenius exponentiation of a discrete 
spectrum of trapping well depths, leading to a discrete hierarchy of time scales. Simi- 
lar effects were found for concentration p = 0.5, while for bias B = 0.5 no oscillations 
were seen. Our result underlines the importance of accounting for intermittent distor- 
tions that will bias the determination of exponents whenever log-periodicity occurs. 
We note that a similar difficulty for determining the fractal dimension of fractal DLA 
clusters in the past has been attributed to the existence of log-periodic modulation 
of the apparent exponent [16]. 

We are indebted to M. Barma and D. Dhar for pointing out an error in a previous 
version and for useful correspondences. We thank the Les Houches Spring School 1997 
"Scaling and Beyond" [17] as well as L. Schafer and his group at Essen University and 
J. Machta and A. Aharony for crucial questions and suggestions. The simulations 
were made at HLR Jiilich in several thousand processor hours. 
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Fig.l: Effective exponent A; = d log / d log t versus time t foi B = 0.98. (combi- 
nation of different runs with different statistics). 

Fig.2: Same as Fig.l for B = 0.999. 

Fig. 3: Same as Fig.l for B = 0.95. 
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